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r  Numerical  simulations  have  been  conducted  using  NRL’s  DIODE2D  computer  code  to  model  the 
steady -state  behavior  of  electron  flow  in  a  radial  diode  and  in  its  adjacent  tri-plate  transmission  line 
(TTL).  Particular  attention  was  paid  to  the  magnitude  of  the  electron  current  flowing  from  the 
cathode  to  the  anode  surface  in  the  TTL.  A  quantitative  value  for  this  effective  power  loss  is  given. 

The  electi^on  current  is  restricted  mainly  to  the  transition  region  in  the  TTL  into  which  there  is  ^seepage" 
of  the  G^|that  is  imposed  in  the  diode  gap.  This  finding  highlights  the  importance  of  that  region  to 
diode  designers.  ^ 
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NUMERICAL  STUDY  OF  ELECTRON  LEAKAGE  POWER  LOSS  IN  A 
TRI-PLATE  TRANSMISSION  LINE 

I.  INTRODUCTION 

There  exist  two  tnaj-'r  alternate  geometries  for  intense  light  ion  beam 
sources  (diodes)  to  be  used  in  light  ion  inertial  confinement  fusion  reactor 
scenarios.  They  are  the  axial  pinch-reflex  diode  (PRD)  pioneered  at  the  Naval 
Research  Laboratory^  and  the  magnetically  insulated  radial  diodes  under 
development  at  Sandia  National  Laboratories^  and  Cornell  University.^  Each  of 
the  two  diode  types  have  their  respective  advantages  and  disadvantages.  The 
low  impedance  (0. 5-2.0  ohm)  PRD's  tested  on  the  Gamble  II  pulsed  power 
generator  at  NRL  proved  themselves  capable  of  delivering  over  70%  of  the  net 
power  flowing  through  the  diodes  to  the  light  ion  (H'*'  and  D'*’)  beams. They 
owed  their  high  ion  production  efficiencies  to  the  theoretically  established 
benefits  of  large  aspect  ratios^  (cathode  radius/A-K  gap)  and  multiple 
reflections  of  electrons  through  the  anode  foil.®  The  geometry  of  the  PRD 
(pinch-reflex  diode)  produced  an  axially  flowing  ion  beam  which  could  then  be 
directed  to  illuminate  a  fusion  pellet  target  located  some  distance  away.^  In 
some  current  fusion  reactor  conceptual  designs  this  can  be  considered  an 
advantage  since  the  ion  beam  source  could  then  be  located  far  enough  away  from 
the  fusion  reactor  vessel  to  protect  them  from  possible  radiation  and  blast 
damage.  In  addition,  there  would  be  no  limit  to  the  number  of  individual 
axial  diode  ion  beam  sources  which  could  be  arrayed  around  a  central  target  to 
deliver  as  much  beam  power  as  necessary  to  the  target.  Unfortunately, 
theoretical  analysis  has  revealed  that  lower  current  ion  beams  lend  themselves 
more  efficiently  to  transport  over  multi -meter  distances.®  To  maintain  a 
given  power  level  while  boosting  transport  efficiency,  higher  impedance  ion 
diode  sources  are  desirable.  Over  the  past  several  years  NRL  has  pursued  an 
experimental  ion  diode  development  program  under  DNA  sponsorship  at  the  Harry 
Diamond  Laboratory  (HDL).  In  that  work,  diode  impedances  range  from  10  to  AO 
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ohms  to  more  closely  match  the  50  ohm  characteristic  impedance  of  the  single 
arm  of  the  AURORA  pulsed  power  generator  used  there  (compared  to  the  1.5  ohms 
of  MRL'S  GAMBLE  II).  To  achieve  the  higher  impedances,  lower  aspect  ratio 
PRO'S  were  employed  and  this  necessarily  lowered  the  diode's  ion  production 
efficiency.  Efficiencies  of  only  20%  to  30%  were  observed.^’  Various 
techniques  are  being  tested  to  improve  this  performance.^^* 

Magnetically  insulated  radial  diodes  offer  a  different  set  of  pros  and 
cons.  They  have  been  shown  to  operate  at  better  than  75%  ion  production 
efficiency.^^  Since  ion  flow  is  focused  to  the  diode  centerline,  beam 
transport  is  not  a  consideration  and  low  impedances  may  be  employed.  This 
also  means  that  target  ignition  will  take  place  inside  the  diode  which  could 
result  in  damage  to  the  device.  The  power  feed  geometry  for  the  radial  diode 
may  also  lead  to  other  difficulties.  Figure  1  presents  a  conceptualization  of 
a  possible  radial  diode  test  apparatus.  As  shown  in  the  cutaway  view,  a  large 
number  of  individually  charged  power  modules  are  arrayed  around  the 
circumference  of  the  device.  They  are  so  designed  to  dump  their  stored 
electrical  energy  simul taneously  into  the  large  disc-like  triplate 
transmission  line  (TTL)  which  fills  the  central  portion  of  the  machine.  This 
triplate  -consists  of  an  anode  disc  sandwiched  between  two  equi spaced  cathode 
discs.  The  interplate  separation  may  also  vary  with  radius.  The  TTL 
(triplate  transmission  line)  terminates  in  the  center  with  an  azimuthally- 
symmetric  radial  diode.  A  possible  configuration  for  that  central  region  is 
illustrated  in  Figure  2.  The  symmetry  of  the  power  flow  to  such  a  diode  is 
one  key  consideration.^^  Another  is  the  highly  inductive  character  of  the  TTL 
itself.  This  high  inductance  can  cause  a  severe  lowering  (via  L  .  dl/dt)  of 
the  voltage  across  the  diode  as  compared  to  the  voltage  output  of  the 
individual  pulsed  power  modules.  Furthermore,  the  resistance  divided  by  the 


Fig.  1  —  Conceptual  reactor  configuration  using  tri-plate 
transmission  line  and  radial  diode  (courtesy,  Sandia  Labora- 


Inductance  of  the  TTL  circuit  can  lead  to  a  phase  shift  between  the  peak 
voltage  across  the  diode  and  the  peak  current  flowing  through  the  diode,  with 
unfortunate  consequences  for  the  power  characteristics  of  the  device.  Those 
problems  have  been  and  are  being  studied  elsewhere.  This  report  focuses  on 
yet  another  possible  TTL  power  loss  mechanism,  namely  that  of  direct  electron 
flow  through  the  vacuum  gap  between  the  plates.  The  properties  of  magnetic 
insulation  are  well  enough  understood  to  ensure  that  an  adequately  insulated 
TTL  could  be  designed  with  sufficiently  large  gap  spacings  for  known  values  of 
diode  current  and  voltage.  However,  there  is  likely  to  be  a  transition  region 
in  the  TTL  in  the  immediate  vicinity  of  the  radial  diode  into  which  will  seep 
a  finite  value  of  the  axial  magnetic  field  imposed  in  the  diode's  A-K  gap. 

The  physics  of  electron  flow  in  such  a  transition  region  cannot  readily  be 
dealt  with  by  analytic  treatment.  For  that  purpose,  this  numerical  simulation 
was  undertaken. 

II.  THE  RADIAL  DIODE  AND  TTL 

The  ideal  diode  and  triplate  transmission  line  geometry  chosen  for 
analysis  is  depicted  in  Figure  3  (drawn  to  scale).  The  cathode  consists  of  a 
solid  shank  5.0  centimeters  in  diameter  and  6.4  centimeters  long  sandwiched 
axially  between  two  infinite  parallel  plates.  The  anode  takes  the  form  of  a 
flat  annular  collar,  mounted  symmetrically  and  coaxially  to  the  cathode 
shank.  The  anode's  inside  diameter  of  5.5  centimeters  yields  a  radial  anode- 
cathode  gap  of  0.5  centimeter  while  its  4.4  centimeter  thickness  allows  a 
clearance  of  1.0  centimeter  from  each  of  the  parallel  cathode  plates.  An 
electrical  potential  difference  of  2.0  megavolts  is  maintained  between  the 
electrodes.  In  addition,  a  field  coil  structure  is  assumed  such  that  a 
constant  and  uniform  background  of  20.0  kilogauss  permeates  the  diode  gap 


4 


V=2.0MV  (Bz)„^^=20kG 

Fig,  3  —  The  radial  diode  and  TTL  modeled 


completely  from  R  =  5.0  to  R  =  5.5  centimeter.  This  field  is  further  assumed 
to  decrease  in  strength  linearly  with  increasing  radius  inside  the  adjacent 
TTL,  until  it  reaches  zero  at  R  =  6.5  cm. 

The  geometry  under  consideration  may  be  subdivided  into  three  separate 
gap  regions  as  shown  in  Figure  4.  Each  half  of  the  "infinite"  triplate  shall 
be  called  Region  I.  It  is  primarily  the  diode  current  flowing  axially  through 
the  cathode  shank  which  generates  an  azimuthal  magnetic  field  in  this  region 
to  insulate  the  electron  flbw.  Region  II  consists  of  the  radial  anode-cathode 
gap  which  is  the  heart  of  the  radial  diode.  There  too  an  insulating  Bg  arises 
from  the  cathode  shank  current.  Note,  however,  that  I^  must  go  to  zero  at  the 
diode's  center  plane.  8.  must  vanish  with  it  there.  For  that  reason,  it  is 

U 

the  strong  applied  axial  magnetic  field,  ,  which  must  be  relied  upon  to 
insulate  the  electrons  in  the  diode.  The  azimuthal  field  component  there 
serves  primarily  to  complicate  the  electron  flow.  Finally,  there  is  the 
transition  volume  of  Region  III.  In  the  lower  third  of  that  region,  §  will 
closely  resemble  that  of  II  but  t  will  be  weaker  and  no  longer  totally 
radial.  In  the  upper  two-thirds  there  is  a  fairly  sharp  B^  gradient  which 
should  act  to  prevent  any  TTL  electrons  from  entering  II.  Exactly  where  those 
reflected  electrons  go  is  the  central  question  treated  herein.  Before 
plunging  into  the  numerics,  it  is  first  illuminating  to  derive  some  zeroeth- 
order  predictions  for  Region  I  and  II  from  simple  analytic  theory. 

Throughout  most  of  Region  II,  the  vacuum  electric  field  has  a  uniform 
value  of  (Ep,  E^)  =  (1.33x10^,  0.0)  statvolts/cm,  while  in  (Ep,  E^)  =  (0.0, 
0.666x10^)  statv/cm.  In  the  absence  of  and  Bg  in  II,  the  expected  space- 
charge  limited  electron  flow  (ignoring  ions)  may  be  approximated  by  the 
Langmui r-Compton^^  expression 
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6  =  in  [  — )  ,  r  =  collector  radius  =  5.5  cm,  Cg  =  emitter  radius  =  5.0  cm, 
e 

L  =  length  of  emitting  shank  =  4.4  cm,  V  =  diode  voltage  =  2.0  iW.  For  the 
given  values  of  radii,  5  =  0.09531  and  =  0.09176.  This  predicts  an 
electron  current  of  3.95  megamperes  in  the  radial  A-K  gap  (N.B.  -  This  ignores 
those  portions  of  the  cathode  shank  in  Region  III).  Significant  electron 
emission  can  also  be  expected  along  the  parallel  cathode  side  plates  of  Region 
I.  The  space-charge  limited  electron  current  density  emitted  there  can  be 
approximated  from  the  relativistic  expression  attributed  to  Friedlander, 
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where  d  =  anode-cathode  plate  separation  =  1.0  cm.  For  the  two  megavolts 
applied,  this  predicts  Jg  =  5.1024  kA/cm^.  In  the  simulation  that  follows, 
electron  emission  is  permitted  on  these  surfaces  between  R  =  5.2  and  R  =  9.5 
cm.  This  yields  an  emission  area  of  198.6  cm^  in  each  TTL  or  a  total  diode 
feed  line  electron  current  of  about  2.0  megampere.  Ignoring  the  effects  of 
magnetic  fields,  this  apparatus  as  modeled  would  therefore  have  a  total 
electron  current  of  about  six  megamperes.  Substituting  an  electron-ion  bi¬ 
polar  curr'  t  in  the  A-K  gap  raises  the  figure  there  by  about  a  factor  of  1.86 
and  thereby  boosts  the  total  diode  current  to  about  9.35  megamperes. 

In  order  to  gauge  the  effect  of  magnetic  insulation  in  Region  I  a 
parallel  can  be  drawn  with  the  analysis  of  the  "critical  current"  in  an 
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idealized  parallel  plate  axial  diode. In  that  device,  the  axial  current 
flow  generates  an  axial  magnetic  field  which  causes  emitted  electrons  to  bend 
toward  the  center  axis.  If  this  B.  (and  therefore  this  I.)  is  large  enough, 
the  electron  will  not  reach  the  anode  plate.  Ignoring  the  effects  of  the 
emitted  electron  space-charge,  a  simple  estimate  can  be  made  for  the  "critical 
current,"  1^.,  necessary  to  produce  an  electron  trajectory  which  just  grazes 
the  anode  surface: 

I^  =  3500  8y  (j)  amperes  (3) 

where  R  =  the  radius  of  emission,  d  =  A-K  gap,  and  $  and  y  have  their  standard 

relativistic  definitions.  Here  8  =  0.98,  y  .  =  4.9,  and  d  =  1 .0  cm  so 

max  max 

that  the  equation  becomes 

I^  =  4.08x10^  R  amperes  (4) 

The  axial  current  taking  part  in  the  physics  of  Eq.  4  may  be  considered  to 
consist  only  of  the  total  cathode  shank  current  (i.e.  -  one  half  of  the  net 
Region  II  current).  By  definition,  the  insulated  electron  flow  in  the  TTL  can 
contribute  no  net  axial  current.  The  previous,  crude  estimate  of  electrons- 
only  Ijj  of  3.95  MA  leads  to  an  axial  current  of  1.98  MA  which  implies,  via 
Eq.  4,  magnetic  insulation  out  to  a  radius  of  about  48.5  centimeters. 

However,  in  a  properly  functioning  radial  diode,  the  current  prediction  of  Eq. 
1  will  not  hold  true.  The  relativistic  electron  gyroradius  in  the  imposed  20 
kG  axial  magnetic  field  is  only  0.41  cm,  well  below  the  0.50  cm  width  of  the 
gap.  Thus,  magnetic  insulation  should  prevail  and  estimates  for  diode 
performance  should  be  sought  directly  from  the  literature  on  magnetically 
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insulated  coaxial  charge  flow.  In  one  of  the  most  recent  comprehensive 
treatments  on  the  subject,^®  it  was  noted  that  the  critical  field  for 
insulation  is  given  oy 


where  r^  is  the  anode  radius,  r^  is  the  cathode  radius 

and  V  =  — ^  (b  .  For  a  =  2  MV,  V  equals  3.9  and  becomes  17  kG,  well  below 

OZ  O  0  0  c 

me 

the  chosen  value  of  B^q  =  20  kG  in  the  diode  under  consideration.  There  also 
appears  in  Reference  18,  rough  parametric  plots  of  ion  current  values  for 
various  combinations  of  diode  characteristics.  From  those  plots,  it  may  be 
crudely  approximated  that  the  ion  current  in  Region  II  should  be  about  2-3 
times  the  Langmui r-Bl odgett  ion  current,  I^^g.  This  I^^g  will  be  equal  to  the 
previously  predicted  3.95  MA  for  electrons  multiplied  by  the  square  root  of 
the  electron-proton  mass  ratio,  mg/m^.  This  leads  to  a  Region  II  ion  current 
of  about  200-300  kA  and  a  half -current  of  100-150  kA  which  in  turn  implies 
magnetic  insulation  in  the  TTL  out  to  a  radius  of  only  2.45  to  3.68  cm.  This 
translates  to  no  insulation  at  all  in  the  TTL  since  its  lower  radius  is  5.0 
cm.  Only  the  full  numerical  simulation  can  accurately  show  what  the  true 
outcome  will  be. 

III.  THE  COMPUTATIOMAL  EXPERIMENT 

The  numerical  simulation  of  the  radial  diode  was  carried  out  with  the 
latest  version  of  NRL's  DI0DE2D  computer  code.^^  This  code  is  two  and  one- 
half  dimensional  (i.e.,  particles  are  allowed  two  spatial  coordinates,  r  and 
z,  as  well  as  three  canonical  momentum  components  Pp,  and  P^  ).  The 
computational  "particles"  are  azimuthal ly-symmetric  charge/current  rings  with 
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charge-to-mass  ratios  corresponding  exactly  to  that  of  electrons  and 
protons.  The  motion  of  these  macroparticles  is  treated  completely 
relativistically.  .0I00E20  is  a  particle-in-cell  (PIC)  code  in  which 
macroparticle  charges  and  current  densities  are  area-weighted  over  the  four 
nearest  grid  points  on  the  computational  mesh.  In  its  electrostatic, 
magnetostatic  treatment,  l)I0DE2D  solves  Poisson's  equation  to 

20 

obtain  (j)(r,z)  and  A.(r,z)  respectively  from  p(r,z)  and  J-(r,z).  From  these 
the  self-consistent  fields  are  given  by  ?  =  -7ij)  and  f  =  (B^,  =  7  x  K. 

The  azimuthal  component  of  the  self-magnetic  field,  B.  ,  is  found  directly 
through  radial  integration  of  Ampere's  Law  over  the  mesh.  Regularly  shaped 
conductor  regions  which  exclude  electric  field  and  magnetic  flux  may  be 
overlayed  on  the  computational  mesh  by  utilizing  a  capacitance  matrix 
technique^^  in  the  Poisson  solving  algorithm.  Tnis  optional  capability  was 
essential  for  this  particular  simulation.  A  static  and  spatially  non-uniform 
axial  magnetic  field  is  imposed  on  the  computational  mesh  in  a  manner 
presented  in  the  treatment  below. 

Figure  5  depicts  the  layout  of  the  computational  mesh  used  to  model  the 
combined  radial  diode  and  TTL  of  Figure  4.  Due  to  its  left-to-right  symmetry, 
only  the  left-half  of  the  device  is  modeled  on  the  computer.  Thus,  the 
current  values  tabulated  in  Section  IV  for  the  numerical  model  must  be 
multiplied  by  two  to  apply  to  the  entire  diode.  In  setting  up  the 
computational  mesh  primary  attention  was  paid  to  ensuring  needed  spatial 
resolution  of  field  and  particle  motions  in  all  regions  of  interest.  As 
presented  in  Section  II,  the  electron  gyroradius  in  the  radial  gap  is 
approximately  0.41  cm.  A  minimum  of  four  cells  per  gyroradius  is  desired 
so  Ar  is  fixed  at  0.1  cm.  Considerations  regarding  the  fast  fourier  transform 
technique  employed  in  the  Poisson  solver  allow  increased  generosity  with  the 
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number  of  axial  gridpoints  so  that  a  az  of  0.05  cm  is  chosen.  For  the 
specified  diode  dimensions  this  yields  a  computational  mesh  spanning  64  cells 
in  z  and  100  cells  in  r. 

The  temporal  resolution  employed  in  the  simulation  plays  an  equally 

significant  role  and,  if  misjudged,  may  lead  to  gross  numerical 

instabilities.  Once  agaiii,  the  electron  gyromotion  is  a  determining  factor. 

Its  frequency,  fg,  is  estimated  to  be  approximately  5.6x10^^  hertz  and  proper 

representation  of  electron  motion  demands  that  f»  .  At  «  1.0.  At  the  same 

time,  it  is  desirable  to  follow  the  simulation  through  several  ion  anode-to- 

cathode  crossing  times  in  order  to  assure  that  a  believable  steady  state  has 

been  reached.  For  the  given  voltage  of  2  MV,  about  1.2x10“^  second  is 

consumed  by  an  ion  traversing  the  radial  gap.  In  order  to  keep  the  cost  of 

the  simulation  at  a  reasonable  level,  this  crossing  time  should  correspond  to 

no  more  than  1000  timesteps.  These  upper  and  lower  brackets  suggest  a 
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timestep  of  At  =  2.0x10  second.  For  this  choice,  f-  At  ~  0.11  and  the  ion 

b  ~ 

crossing  time  equals  600  timesteps. 

Particles  are  injected  into  the  diode  gap  using  the  space  charge  limited 
emission  condition.  (Just  enough  charge  is  emitted  so  that  the  component  of 
the  electric  field  normal  to  the  emission  surface  vanishes.)  In  this 
simulation,  electron  emission  is  permitted  up  to  a  radius  of  9.5  cm  as  well  as 
all  along  the  axial  cathode  surface  to  within  0.2  cm  of  the  diode  center 
plane.  Emission  was  suppressed  at  the  extrema  of  the  computational  mesh 
because  of  expected  distortion  of  field  values  near  the  boundaries.  In 
contrast  ion  emission  is  restricted  to  the  inner  axial  anode  surface  except 
for  the  region  with  0.2  cm  of  the  diode  center  plane.  Emitting  sufaces  are 
indicated  by  the  heavy  lines  in  Figure  5.  Since  only  the  left  half  of  the 


diode  is  being  simulated  complete  symmetry  demands  that  any  particles 
impacting  the  right-hand  computational  boundary  are  perfectly  reflected  back 
into  the  simulation  mesh. 

The  question  of  electric  and  magnetic  boundary  conditions  is  not  trivial 
to  handle  properly.  Since  no  particle  emission  is  permitted  closer  than  0.5 
cm  (5  data  cells)  to  the  outer  radial  boundary  of  the  TTL,  there  is  no  serious 
problem  in  assuming  no  change  to  the  source-free,  static  potential  solution 
there  during  the  entire  simulation.  For  that  reason,  the  electrostatic 
potential,  rises  linearly  in  z  from  o  to  across  the  1  cm  gap  along  that 
boundary  (see  Figure  4).  (The  magnetostatic  boundary  condition  there  is 
irrelevant  since  there  is  no  imposed  in  Region  I.)  The  right  boundary  of 
Region  II  is  more  difficult.  It  clearly  does  not  remain  source-free  during 
the  simulation.  Unfortunately,  the  capacitance  matrix  in  the  Poisson-solver 
must  be  recalculated  every  time  the  boundary  values  of  the  potentials 
change.  This  is  a  costly  procedure.  Rather  than  pay  the  price  of  matrix 
recalculation,  it  was  decided  to  freeze  the  values  of  ij)  along  that  boundary  as 
well  at  their  source-free,  logarithmically-spaced  values  (as  in  Figure  4).  In 
addition,  no  self  magnetic  flux  penetration  across  that  diode  centerplane  is 
permitted.  These  conditions  do  some  violence  locally  to  the  f  and  ^  field 
structure  there  but  the  gross  dynamics  of  Region  II  should  remain  unaltered. 

IV.  RESULTS 

To  better  understand  the  physics  of  the  device,  two  distinct  simulations 
were  conducted.  In  the  first,  a  benchmark  for  the  radial  diode  operation  was 
sought  by  permitting  electron  emission  only  on  the  cathode  shank.  This 
provided  better  estimates  for  the  ideal  magnetic  insulation  to  be  expected  in 
the  TTL  region  and  it  gave  a  reference  figure  for  the  maximum  operating  power 


of  the  diode  itself.  The  second  simulation  allowed  electron  emission  all 
along  the  cathode  surface  in  the  TTL  as  well  as  along  the  shank.  Sample 
particle  position  plots  were  generated  to  depict  the  nature  of  the  electron 
flow  there  in  conjunction  with  plots  of  the  radial  profile  of  electron  current 
densities  striking  the  anode  surface  in  the  TTL.  Direct  measurement  was  made 
of  the  electron  current  flowing  in  Region  I  to  yield  the  percentage  of  total 
power  which  never  reached  the  diode. 

Case  I,  which  simulated  only  the  radial  diode,  was  run  over  1200 
timesteps  of  2.0x10“^^  second  each.  The  electron  and  ion  currents  of  the 
half -diode  being  modeled  are  depicted  in  Figure  6.  (Multiply  by  2  to  obtain 
full  diode  currents.)  During  the  first  400  timesteps,  only  electron  (no  ions) 
were  permitted  in  the  system.  The  dots  on  the  plot  represent  emitted  electron 
current.  The  solid  lines  are  the  collected  ion  and  electron  currents  as 
indicated.  Note  that  the  electron  flow  was  completely  insulated  in  the  diode 
during  electron-only  operation.  A  steady  electron  current  of  about  75  kA 
continued  to  be  emitted  by  the  cathode  shank  but  an  equal  electron  current  was 
continuously  absorbed  by  the  shank.  This  is  characteristic  of  ideal  magnetic 
insulation. A  constantly  replenished  cloud  of  electrons  is  formed 
adjacent  to  the  cathode  surface  and  extends  some  finite  radial  distance  into 
the  A-K  gap.  This  is  illustrated  in  the  sample  electron  position  plot  of 
Figure  7.  (Note  again  that  only  half  of  the  diode  is  shown.  The  other  half 
is  completely  symmetric  about  the  center-plane.)  Mo  net  current  is  flowing 
through  the  diode.  Therefore,  the  TTL  is  completely  uninsulated.  Electrons 
in  that  region  are  free  to  migrate  toward  the  anode  limited  only  by  the  self- 
fields  they  may  generate. 

At  T  =  400  At  ,  ion  emission  was  "turned  on"  all  along  the  cylindrical 
anode  surface.  The  resultant  neutralization  of  electron  space  charge  caused  a 
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sharp  increase  of  electron  emission  as  seen  in  Figure  6.  In  the  new 
equilibrium  state  formed  by  the  two-specie  flow,  almost  3  MA  was  continuously 
emitted  and  reabsorbed  by  the  cathode  half-shank  to  maintain  the  enlarged 
negative  space-charge  cloud  depicted  in  Figure  8.  The  densest  portion  of  the 
cloud  is  now  within  one  millimeter  of  the  anode  surface.  This  accounts  for 
the  steady-state  ion  half-current  of  about  1.1  MA,  far  above  the  simple 
Langmui r-Bl odgett  prediction  of  100-150  kA.  In  contrast,  the  net  electron 
half-current  is  a  mere  250  kA  compared  to  the  nearly  2  MA  L.-B.  estimate.  The 
net  current  flowing  through  one-half  of  the  diode  is  thus  about  1.35  MA. 
According  to  Eq.  4,  this  should  provide  magnetic  insulation  in  the  TTL  out  to 
a  radius  of  about  33  cm  if  the  imposed  axial  magnetic  field  were  not 
present.  However,  B^  cannot  be  ignored.  It  is  the  key  to  ion-efficient 
operation  of  the  radial  diode.^^  In  combination  with  tHe  diamagnetic  self- 
fields  generated  by  charge  circulation  in  the  electron  cloud,  the  imposed  B^, 
forms  an  axial  magnetic  field  in  the  A-K  gap  whose  '’verage  radial  profile  is 
as  shown  in  Figure  9.  The  net  field  gradient  is  jfti  u  "H  the  electrons-only 
case  but  exceeds  140  kG/cm  when  full  two-specie  flow  is  present.  The  field 
strength  peak  of  almost  70  kG  near  the  anode  surface  in  Region  II  acts  as  an 
effective  barrier  against  electron  penetration  there.  The  details  of  the  role 
this  B-field  plays  in  the  A-K  gap  physics  is  discussed  elsewhere. Of 
primary  importance  here  is  the  "spill-over"  of  this  field  into  Region  II,  the 
lower  end  of  the  TTL. 

In  the  model  chosen  for  this  simulation,  the  imposed  20  kG  field  is 
specially  uniform  between  R=5.0  and  R=5.5  cm.  It  decays  linearly  to  zero  with 
increasing  radius  between  R=5.0  and  R=6.0  cm  in  the  TTL.  To  that  background 
field  must  be  added  the  particle  self-fields.  Sample  values  of  the  net 
B2(r,z)  are  presented  in  Figure  10  for  the  cases  with  and  without  ions  .  Note 
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Fig.  8  —  Sample  ion  and  electron  positions  for  the  full, 
equilibrium  “Benchmark”  case 
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Fig.  9  —  Radial  profile  of  in  the  A-K  gap  (averaged  in  the 
z  *  direction)  for  the  “Benchmark"  case 
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Fig.  10  —  Sample  values  of  near  the  anode  for  the 
pure  nKlial  diode  “Benchmark”  case 
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that  the  Region  II  fields  are  hardly  modified  at  all  in  the  electrons-only 
steady  state.  The  addition  of  ions  into  the  system,  as  already  pointed  out, 
does  make  a  significant  difference.  Particularly  on  the  anode  side  of  Region 
III,  the  values  of  are  elevated,  resulting  in  a  somewhat  steeper  magnetic 
gradient.  The  modification  of  the  fields  there,  though  significant,  is  quite 
small  in  comparison  to  the  tripling  of  inside  Region  II.  Therefore,  it  is 

the  imposed  field  gradient  in  the  TTL  that  will  be  responsible  for  loss  of 
insulation  there.  Over  most  of  the  axial  distance,  that  gradient  is  little 
changed  by  current  flows  in  the  diode. 

The  second  simulation  run.  Case  II,  tested  the  full  combination  of  the 
radial  diode  and  the  TTL.  Figure  11  provides  the  corresponding  point  plot  of 
emitted  electron  current  as  well  as  traces  of  collected  ion  and  electron 
currents.  Once  again,  the  numerical  model  was  run  for  the  first  400  timesteps 
without  any  ions  present.  The  plot  shows  about  620  kA  of  electron  current 
continuously  emitted  in  the  steady  state  with  a  full  450  kA  of  it  being 
collected  by  the  anode.  This  leaves  about  170  kA  of  electrons  to  sustain  any 
recirculating  electron  cloud  in  Region  II  or  in  the  TTL  (compared  with  75  kA 
in  Case  I).  The  fate  of  the  450  kA  net  flow  can  be  deduced  from  the  sample 
electron  position  plot  of  Figure  12.  Electrons  appear  to  be  bridging  the  TTL 
gap  out  to  a  radius  of  about  7.5  centimeters  or  more.  The  fact  that  electrons 
make  it  across  at  all  is  consistent  with  the  previously  determined  lack  of 
electron-only  current  flow  across  the  radial  diode  gap.  .Any  B.  insulation  in 
the  TTL  must  therefore  arise  from  the  electron  current  bridging  the  gap 
there.  A  simple  estimate  may  be  obtained  for  the  minimum  insulation  radius  by 
equating  the  critical  current  of  Eq.  4  with  the  net  TTL  current  below  a 
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Fig.  11  —  Electron  and  ion  currents  in  the  full  radial  diode-TTL  case 


Fig.  12  —  Sample  electron  positions  with  electrons-only  flow  in 
the  diode-TTL  case 
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radius,  R,  assuming  a  uniform  Jg  of  5.1  kA/cm^  as  determined  from  Eq.  2.  This 
produces  the  expression: 

TT  (r^  -  5.5^)  5.1x10^  =  4.08x10^  R  (6) 

which  yields  a  minimum  radius  of  about  6.9  centimeters,  in  agreement  with  what 
is  observed  in  the  simulation.  For  a  more  complete  picture  of  the  electron 
flow  pattern  in  the  TTL,  the  emitted  and  collected  electron  current  densities, 
Jg,  were  plotted  in  Figure  13.  An  integration  of  the  collected  current  yields 
425  kA,  which  accounts  for  essentially  all  of  the  net  450  kA  of  electron 
current  flowing  through  the  device.  Furthermore,  the  primary  peak  of 
collected  Jg  lies  below  6.5  cm  in  radius  in  agreement  again  with  the 
theoretical  estimate  of  6.9  cm.  It  is  also  interesting  to  note  the  profile  of 
the  emitted  jg.  As  expected,  the  highest  peak  lies  above  R=9.0  cm  where  the 
electron  emission  region  begins  and  the  surface  E-fields  are  highest  due  to 
two-dimensional  distortions  in  the  equipotential  line  there.  At  that  high 
radius,  however,  the  8.  is  more  than  adequate  to  prevent  tiiu  emitted  electrons 
from  reaching  the  anode.  They  are  bent  around  to  strike  the  cathode  again  at 
around  R=8.0  cm,  causing  a  virtual  shut-off  of  emission  there  due  to  a  buildup 
of  negative  space-charge.  The  emission  profile  is  much  more  steady  below 
R=7.5  cm  where  B„  is  not  strong  enough  to  confine  electron  flow  near  the 
cathode. 

As  in  the  previous  case  which  looked  only  at  the  radial  diode,  ion 
emission  in  Region  II  was  "turned  on"  at  T=400  At.  The  relaxation  of  the 
system  to  a  new  steady  state  can  be  traced  in  Figure  11.  By  T=1200  At  ,  the 
ion  half-current  collected  has  settled  to  about  925  kA  while  that  of  the 
electrons  is  about  300  kA.  Figure  14  depicts  a  sample  set  of  electron  and  ion 
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positions  at  equilibrium.  The  electron  plot  is  very  different  from  that  of 
Figure  12.  As  previously  stated,  an  ion  half-current  of  about  1.0  MA  is  more 
than  adequate  to  insure  magnetic  insulation  throughout  the  TTL  in  the  absence 
of  any  B^.  Thus  it  is  no  surprise  to  see  electrons  there  confined  to  a  thin 
sheath  near  the  cathode  surface  down  to  about  6.5  cm,  where  the  applied  axial 
B-field  begins  to  rise  in  value  from  zero  to  a  full  20  kG  at  R=5.5  cm.  In 
that  annular  band  between  5.5  and  6.5  cm,  a  very  significant  leakage  of 
electron  current  to  the  anode  occurs.  This  shorting  should  not  manifest 
itself  were  it  not  for  the  applied  B^  gradient  there.  To  quantify  the  amount 
of  power  lost  due  to  this  leakage,  the  radial  profiles  of  the  electron  current 
density  emitted  and  collected  in  the  TTL  are  plotted  in  Figure  15.  The 
emission  profile  is  quite  irregular.  Peaks  and  troughs  alternate  down  the 
entire  length.  Enhanced  emission  at  one  point  leads  to  negative  charge 
buildup  downstream  in  the  electron  sheath,  causing  reduced  emission  there.  On 
the  other  hand,  the  collected  Jg  plot  consists  of  a  single,  intense  band  below 
R=6.5  cm  as  could  be  guessed  from  the  sample  position  plot  of  Figure  14.  A 
total  of  about  205  kA  of  electron  current  is  impacting  there  in  the  TTL.  This 
corresponds  to  a  full  two-thirds  of  the  net  electron  current  flowing  through 
the  entire  device.  What  is  even  more  disturbing  is  that  0.41  terawatts  out  of 
a  total  of  2.45  TW  flowing  through  the  half -diode  modeled  are  being  lost  in 


V.  CONCLUSIONS 


The  major  findings  of  this  set  of  numerical  simulations  are  summarized  in 
Table  1.  As  stated  in  the  introduction,  the  optimization  of  net  ion 
production  efficiency  is  a  central  objective  of  pulsed  power  diode  research. 
The  table  shows  a  loss  of  almost  20%  of  the  net  input  power  in  the  TTL  for  the 
configuration  modeled.  The  loss  can  be  attributed  to  the  steep  gradient  in 
in  the  transition  region  between  the  TTL  and  the  radial  diode.  The  result  is 
about  an  8%  drop  in  the  total  ion  production  efficiency  of  the  TTL-diode 
combination.  That  can  probably  be  tolerated  for  most  practical 
applications.  However,  any  experimental  configuration  different  from  that 
shown  in  Figure  3  will  experience  a  different  and  possibly  greater  power  loss 
in  the  TTL.  One  must  clearly  exercise  care  in  setting  the  TTL  gap  and  in 
shaping  any  imposed  magnetic  field  profiles  in  the  transition  region. 


Table  1  —  Sununary  of  results 


CASE 

I 

CASE 

II 

ions  present  ? 

no 

yes 

no 

yes 

^^e^emitted 

0.15  MA 

6.0  MA 

1.24 

MA 

5.0  MA 

^^e^collected  in  TTL 

- 

- 

1 

0.9 

MA 

0.4  MA 

^ ^e^collected  in  diode 

- 

0.5  MA 

0 

0.2  MA 

^^1^ collected  ^ 

1 

- 

2.2  MA 

- 

1.85  MA 

- 

0.815 

- 

0.902 

P 

Ld  p„  + 

D  TTL 

1.0 

■ 

0.833 

REFERENCES 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 


12 


13 

14 

15 

16 

17 

18 
19 


S.  J.  Stephanakis,  D.  Mosher,  S.  A.  Goldstein,  et.al..  Bull.  Am.  Phys. 
Soc.  1130  (1977). 

O.  J.  Johnson,  G.  W.  Kuswa,  A.  V.  Farnsworth,  Jr.,  J.  P.  Quintenz, 
et.al.,  Phys.  Rev.  Lett.  610  (1979). 

P.  Dreike,  C.  Eichenberger,  S.  Humphries,  and  R.  Sudan,  J.  Appl.  Phys. 

85  (1976). 

S.  J.  Stephanakis,  J.  R.  Boiler,  G.  Cooperstein,  S.  A.  Goldstein, 
et.al..  Bull.  Am.  Phys.  Soc.  _23,  907  (1978). 

S.  A.  Goldstein  and  R.  E.  Lee,  Phys.  Rev.  Lett.  _35,  1079  (1975). 

J.  W.  Shearer,  Lawrence  Livermore  Lab  Report  UCRL-52129  (1976). 

P.  F.  Ottinger,  S.  A.  Goldstein,  and  0.  Mosher,  NRL  Memo  Report  No.  4548 
(1981). 

0.  G.  Colombant,  S.  A.  Goldstein,  and  D.  Mosher,  Phys.  Rev.  Lett.  45, 
1253  (1980). 

S.  J.  Stephanakis,  J.  R.  Boiler,  G.  Cooperstein,  et.  a1.,  Bull.  Am. 

Phys.  Soc.  1031  (1979). 

R.  A.  Meger,  F.  C.  Voung,  A.  T.  Orobot,  et.  a1.,  NRL  Memo  Report  4477 
(1981). 

R.  J.  Barker,  S.  A.  Goldstein,  and  A.  T.  Orobot,  NRL  Memo  Report  No. 

4642  (1981). 

R.  J.  Barker  and  S.  A.  Goldstein,  to  appear  as  a  NRL  Memo  Report  (1982). 

Particle  Beam  Fusion  Progress  Report  April  1978  through  December  1978, 
Sandia  National  Labs,  SAND  79-1011  (1  979). 

D.  8.  Seidel  and  J.  P.  VanDevender  in  Particle  Beam  Fusion  Progress 
Report  June  1979  -  December  1979,  Sandia  National  Labs,  SAnd  80-09)4 
T19807: -  - 

I.  Langmuir  and  K.  T.  Compton,  Rev.  of  Modern  Physics  _3,  191  (1931). 

F.  Friedlander,  R.  Hechtel,  H.  Jory,  et.al.,  Megavolt-Megampere  Electron 
Gun  Study,  DASA  2173,  Varian  Associates,  Palo  Alto,  tA  (1968). 

Robert  K.  Parker,  Technical  Report  AFWL-TR-73-92,  U.  S.  Air  Force 
Weapons  Laboratory,  Albuquerque,  NM  (1  973). 

K.  D.  Bergeron,  Phys.  Fluids  693  (1977). 

R.  J.  Barker,  A.  T.  Orobot,  R.  E.  Lee,  and  S.  A.  Goldstein,  Proc.  9th 
Conf.  on  Numerical  Simulation  of  Plasmas,  Evanston,  IL  (1980). 


31 


20 


R.  J.  Barker,  Banach  Center  Publication's,  255;  Warsaw,  Poland  (1  975). 

B.  L.  Buzbee,  F.  W.  Dorr,  and  J.  A,  George,  Technical  Report  CS-71-195, 
Computer  Science  Dept.,  Stanford  University,  Stanford,  CA  (1971). 

R.  N.  Sudan  and  R.  V.  Lovelace,  Phys.  Rev.  Lett.  117A  (1973). 

R.  V.  Lovelace  and  E.  Ott,  Phys.  Fluids  17,  1263  (1974). 

A.  Ron,  A.  A.  Mondelli,  and  N.  Rostoker,  IEEE  Trans.  Plasma  Science,  PS- 
U  85  (1973). 

R.  J.  Barker,  S.  A.  Goldstein,  and  A.  T.  Drobot,  IEEE  Int.  Conf.  on 
Plasma  Science,  Madison,  WI  (1980). 


32 


